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FLAT-PLATE COLLECTOR PERFORMANCE EVALUATION 
THE CASE FOR A SOLAR SIMULATION APPROACH 

By 

F. F. Simon and Paul Harlamert 


ABSTRACT 

A method is proposed for determining the performance of 
flat-plate solar collectors using a simulated t! sun" . Collector 
test variables that will help establish the basis for the in- 
door test facility at the Lewis Research Center are discussed. 
The use of the indoor testing should permit a standard test for 
the convenient and accurate determination of collector perform- 
ance. Preliminary test results are reported as an example of 
the type of collector performance data to be expected from the 
simulation approach. 



FLAT-PLATE COLLECTOR PERFORMANCE EVALUATION 


THE CASE FOR A SOLAR SIMULATION APPROACH 

By 

F. F. Simon and Paul Harlamert 
Lewis Research Center, Cleveland, Ohio 44135 


INTRODUCTION 

The NASA LeRC is investigating the possible application of 
solar energy as a clean, nondepleting source to help our nation’s 
energy needs. This investigation includes the generation of 
electricity, the heating and cooling of buildings and the produc- 
tion of clean fuel (ref. .1) . Approximately 25 percent of the 
nation’s total energy is used for heating and cooling of homes 
(ref. 2) . Solar energy captured by flat-plate collectors have 
been and are presently being used to heat homes. To make a 
significant impact on our energy needs in this area, efficient, 
economical, reliable solar collectors must be developed. 

The experimental determination of collector efficiency is 
an important requirement in the evaluation of a collector de- 
sign. Much of the past collector performance data are of 
limited value because the data are applicable only for a given 
location under a given weather condition. This limitation does 
not permit a standardized comparison between various collector 
designs. One approach to get around this difficulty is the simul- 
taneous testing of different collector designs (ref. 3) . This 
method, however, has practical limitations. Another approach is 
to use the "standard test" suggested by Whillier (ref. 4) in which 
the data collected are used to calculate performance parameters. 

This approach was used by Gupta and Garg (ref. 5) for the compari- 
son of various collector designs. Because of the variation of 
such parameters as wind speed, air temperature, etc. , the perform- 
ance data of Gupta and Garg show much scatter: What is needed for 

an evaluation and comparison of collector performance is a stand- 
ard method under controlled environmental conditions. This paper 
attempts to demonstrate that the requirements for standard 
collector performance tests and the determination of key perform- 
ance parameters are best satisfied by a solar simulation approach. 
After a discussion of the collector testing requirements a descrip- 
tion is given of the facility presently being built up at the Lewis 
Research Center for the indoor testing of collectors under simulated 
conditions. Preliminary tests on a prototype facility system are 
reported to help validate the solar simulation approach. 
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COLLECTOR TESTING REQUIREMENTS 

A discussion of testing requirements is a useful exercise 
in helping to determine the variables that require controls. In 
deciding on the testing approach, two goals of the testing pro- 
gram need to be kept in mind. It is necessary to have standard 
conditions that will permit a comparison of different collectors, 
and it is required that collector performance parameters be 
determined that will allow a prediction of collector performance 
under other than standard conditions. This discussion will use a 
collector model of the traditional design shown in figure 1. This 
collector design has channels for the heat exchange medium that 
collects the energy absorbed on the plate, one or more windows to 
allow the entry of solar radiation and reduce convective and radia 
tion losses and insulation to reduce rear and side conduction 
losses . 

Collector Performance Variables 


The basic equation relating collector performance to the 
absorber plate temperature is as follows: 


1 = <W e - V T p - T a ) /q T « 

Calculation of the collector efficiency ( r) ) by equation 1 re- 
quires a knowledge of the product of absorptivity ( OC ) and trans- 
mittance ( T e ) , the overall heat loss coefficient ( Ql ) , the 
plate and ambient temperatures ( Tp,T a ) and the amount of solar 
radiation ( qx ) • In the practical determination of collector per- 
formance, it is more convenient to measure the temperature of the 
fluid as it enters and leaves the collector rather than measuring 
the plate temperature. Therefore, equation 1 is neither practical 
or convenient for the determination of collector efficiency. How- 
ever, as is shown later, if an experimental determination of the 
product of absorptivity and transmittance ( 0;'T e ) is required, 
then equation 1 is needed. Whillier (ref. 6) expressed the collec- 
tor efficiency in terms of the average fluid temperature and a 
plate efficiency factor ( F T ) as follows: 


n = Ft ' [«T e - U L (T f - T a )/q x ] (2) 

An even more convenient expression of collector efficiency is given 

in reference 6 in terms of the inlet fluid temperature and a factor 

called the heat removal efficiency factor (F d ), as follows: 

k 

" - F R L m e - “l < T 1 - T P /q l] (3) 
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The heat removal efficiency factor is a function of flow 
rate (G) , heat loss coefficient (Ul ) , plate efficiency factor 
(F f ) and fluid heat capacity ( C p ) . 

F r = f (G, U L , F», C p ) C4) 

For a given collector design, the plate efficiency factor is ex- 
pressed in terms of fluid convective coefficient ( h p ) and the 
heat loss coefficient. 1 


F* = f(h f , U L ) CS) 

The product of absorptivity and transmittance is a function of 
the angle of incidence ( 6^ ) and the relative quantities of 
diffuse ( q ) and direct radiation ( q dr ) . 


aT e = f(e i’ q df’ c ‘dr ) C6 

For a given collector, the overall heat loss coefficient is a 
function of the absorber plate temperature (T p ) , the ambient tem- 
perature (T a ) and the exterior convective coefficient or what 
is commonly called the wind coefficient ( h Q ) . 


U L - fCT p . T a , h o ) (7) 

For an accurate evaluation of collector heat loss, the ambient 
temperature requires a humidity corrective to determine the "effec- 
tive" sky temperature (T s ). The effective sky temperature is 
different for a clear day than for a cloudy day. Equation 7 is 
more correctly stated as: 


Ul = f(T p , T s , h 0 , T a ) (7a) 

The amount and spectral quality of solar radiation is a func- 
tion of such things as air mass, cloud cover, air pollution, water 
vapor, dust and the molecular components of the atmosphere. All 
the above functional relationships indicate the large number of 
variables which can affect the performance of a given collector 
design. In summary, the collection efficiency of a given collec- 
tor design is a function of the following: 
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1. fluid flow rate 

2. ambient temperature 

3. plate temperature 

4. inlet fluid temperature 

5. wind speed 

6. amount ancl spectral quality of solar radiation 

7. incident angle of solar radiation 

8. relative quantities of diffuse and direct radiation 

Consideration of the above variables gives the basis for a 
standard testing approach. For a method that can be used for 
comparing different collector designs, we need to maintain con- 
stant values of the flow rate, ambient temperature and wind speed. 
With these standard conditions we can vary the inlet fluid tem- 
perature and the amount of thermal radiation so as to generate 
standard performance curves. Such performance curves allows us 
to compare collector performance for a given inlet temperature 
and value of thermal radiation. The standard method should em- 
ploy a thermal radiation source which is equivalent to the spec- 
tral output of air mass 2. Air mass is approximately equal to 
sec 0 Z , where & z is the zenith angle of the Sun. Therefore, 
a normal path of the Sun through the atmosphere is air mass 1. 
Benning (ref. 7) has suggested that the "average” sunlight for 
the North American continent is best represented by air mass 2 
spectra. The standard radiation source should give a fixed ratio 
of direct to diffuse radiation (or all direct) and permit a con- 
trol of the incident angle of radiation. A universal agreement on 
the above considerations for standard testing conditions would 
allow, in addition to comparative testing at one site, a basis for 
the comparison of collector test results made in different parts 
of the world. In addition, since th e solar energy work at the LeRC 
will also involve long-range outdoor testing of collectors, it is 
important to have a testing standard that can be used for judging 
the performance life of a collector. This can be determined from 
indoor tests prior to and after the outdoor tests. 

Collector Performance Parameters for Standard Tests 

With the standard test method described above and a well in- 
strumented collector, it is possible to experimentally determine 
the collector performance parameters of equations 1, 2 and 3. 

From the slopes and intercepts of the plots of collection efficiency 
versus the ratio of temperature difference to heat flux, the follow- 
ing parameters may be determined: 
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These parameters are "indices of performance" of a collector 
and except for the heat loss coefficient ( Ul ) should be made as 
high as possible. These indices become important in the comparison 
of different collector designs as discussed by Whillier (ref. 6) . 
Since the collector plate efficiency factor ( F f ) is a measure of 
how collector geometry effects performance, we can use this factor 
for comparing designs if pumping costs are not a consideration. 
Whillier (ref. 6) states that the best collector design is one 
which has the lowest value of the quotient dollars/(square feet-P ) . 
Because the heat-removal efficiency factor is related to flow rate, 
as well as to collector geometry effects, it can be used for com- 
paring collector designs when pumping costs are a consideration. 
Whillier (ref. 6) states if flow rates are used which give the 
same pumping costs for all collectors then the best collector de- 
sign is the one that has the lowest value of the quotient dollars/ 
(square feet-FR ) . In addition to the collector design considera- 
tions of Whillier, if a collector is to perform well at high tem- 
peratures, then it is necessary that it have low heat loss (IowUt ) 
and as high a value of transmittance ( T e ) and absorptivity ( a ) 
as possible. Since these "indices of performance" will be deter- 
mined for the standard conditions mentioned above, it should be 
possible to conveniently select the best collector designs for a 
given system. 

Another important collector performance parameter is the 
collector heat capacity which is used for evaluating collector 
warm-up time and the effect of clouds. The heat capacity is 
easily determined with a solar simulator and a well instrumented 
collector. The effect of cloud cover may be investigated by vary- 
ing the power to the simulator so as to simulate the variation of 
the radiant flux caused by clouds. 

Collector Performance Parameters for Outdoor Tests 


The two efficiency factors (F'&F^) determined from the 
standard tests are, strictly speaking, valid only for the condi- 
tions of the standard tests. However, since they Is ve a small 
dependence on such variables as temperature, solar radiation and 
wind velocity, the parameters can be used for the prediction of 
collector efficiency for a wide range of environmental conditions. 
Since the collector heat loss coefficient is dependent on the 
ambient temperature (approximately 1 percent change for every 10 
degree change in T a ) , the question arises as to what heat loss 
coefficient can be expected for ambient temperatures far different 
than that used in the standard tests? A possible way to answer 
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this is to separate out in the standard tests the heat loss due 
to radiation and convection from the overall heat loss. With a 
knowledge of the radiation properties of the collector it be- 
comes possible to calculate the radiation losses and correlate 
the convection losses. If we take a simple single window collec- 
tor as an example, where the heat loss can be expressed as follows: 


where 


h (T - T ) + ^ ~ 

C P § 

€ £ 

P g 


oCT' 4 - t' 4 ) 

v p g 


( 8 ) 


h = c (T - T f 
c P g 


(8a) 


tn in 

q T = h n (T - T ) + e a (T - T ) 
L ° v g a y g v g 3. J 


( 9 ) 


To solve for the heat loss for a given wind condition, ambient 
temperature and plate temperature, we need to solve equations 8 and 
9 for three unknowns ( q^ , T„ , h c ) . By using a well instrumented 
and well insulated collector in the standard tests, we can determine 
the total heat loss of equations 8 and 9; and, after making a calcu- 
lation for the radiation loss of these equations, a correlation can 
be made of the convective heat transfer. This correlation of the 
convective heat transfer coefficient (equation 8a) will yield the 
constants c and n for a given collector design. With this informa- 
tion the calculation approach of equations 8 and 9 can be used for 
determining the heat loss and heat loss coefficient ( Uj-^q^/Tp-Tg ) 
at other ambient temperatures. 

If the standard tests utilize a thermal radiation source with 
all direct radiation then the value of the product of absorptivity 
times transmittance ( Org) obtained from the standard tests is 
valid for the portion of solar radiation which is direct. The 
value of this product ( cti-Cg ) to be used for diffuse radiation may 
be obtained by integrating with respect to the radiation incident 
angle ( ©i ) the value of the product 0!T e determined indoors 
for different angles of incident. 


INDOOR FACILITY FOR COLLECTOR TESTING 

The indoor test facility at the Lewis Research Center is so 
designed that the test requirements outlined in the previous sec- 
tion can be met. Figure 2 gives a view of the facility with the 
collector in a position to receive radiation from the simulator 
at a zero incidence angle. The simulator is designed (see ref. 8 
for details). for air mass 2 type solar radiation. It consists of 
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143 tungsten-Halogen 300 watt lamps (43 kW) placed in a modular 
array with fresnel lenses placed at the focal distance so as 
to collimate the radiation. A blower at the top of the facility 
prevents overheating of the lamps. Because the lamps have a 
dichroic coating to reduce the infrared output, it is possible to 
get a good spectral simulation of the Sun f s radiation. A calcu- 
lation of the absorptivity of a selective surface under simulated 
conditions indicates an error of only 0.1 percent in cc when 
compared to air mass 2 radiation. There is virtually no change in 
the spectral output as a result of varying the power of the lamps 
which can deliver radiant fluxes up to 320 BTU/hr-ft. 

The simulator which is normally fifteen feet from the collector 
provides a collimated beam that covers an area four feet by four 
feet with a flux uniformity of +5%. It was determined that a great 
economic advantage could be achieved if some latitude were permitted 
on the degree of collimation of the simulator. The collimation half 
angle is 4^ degrees, which is judged to be adequate for flat-plate 
collector testing. For normal incident tests, the error due to im- 
perfect collimation is small and does not become noticeable for angles 
of incident up to sixty degrees. Beyond an incident angle of sixty 
degrees, flat-plate collectors are not effective in collecting 
solar energy and therefore testing beyond this incident angle is 
not required. Controls have been provided which permit an adjust- 
ment of the radiation plane of the simulator so as to accommodate 
different collector tilt angles. 

The collector to be tested is placed on a stand (fig. 2) 
which allows for adjustment of the collector tilt angle and varia- 
tion of the incident angle of radiation. The base of the collector 
stand is on a turntable that rotates to simulate the rotation of 
the Barth with respect to the Sun. This allows for variation of 
the angle of incidence. The angle of incidence can also be varied 
by changing the tilt angle of the collector. This procedure has 
the disadvantage that the collector convective heat loss varies 
slightly with the angle of tilt and thus adds another performance 
factor for comparing different collector designs. The purpose of 
the fan shown blowing on the collector in figure 2 is to provide 
controlled convective thermal losses that would be created by a 
wind of constant speed. At one end of the facility large doors 
are available which can be opened for outdoor heat loss determina- 
tions . 

For the testing of collectors utilizing water systems, a 
thermal storage unit is provided. The water contents can be 
flowed through a heat exchanger (fig. 3) for quick adjustment of 
the inlet temperature to the collector. Without the heat exchanger, 
it would be required to vary the temperature of the water in storage. 
This would involve long waiting times. In collectors utilizing air 
as a heat exchange medium, the air is sent through a heat exchanger 
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to add thermal energy for control of the inlet temperature (fig. 4-) . 
In addition to evaluating collector efficiencies, fluid pressure 
drop within the collector is also. evaluated to determine pumping 
costs. 


Data from a well instrumented collector including glass tempera- 
tures, plate temperatures, inlet and outlet fluid temperatures, flow 
rates, pressure drop, etc., are recorded on tape with an automatic 
voltage digitizer system. This is done because of the large number 
of tests and measurements to be made. The data can be quickly and 
accurately reduced with a high-speed computer. Some of the key data 
will also be recorded with pen recorders. The thermal radiation 
flux will be measured by using a water-cooled pyrheliometer. 


PRELIMINARY COLLECTOR TESTS 


A Apparatus 

Preliminary tests were run with a small version (prototype) 
of the solar simulator to check out the test equipment and the 
test approach. The collector used in the test was one foot by one 
foot insulated by styrofoam and with one cover glass (fig. 5) . 

The absorber plate was made by soldering 3/8" o.d. copper tubes to 
a 0.03-inch copper plate with the copper tubes being connected to 
two 3/4— inch copper headers which were used for flow distribution. 

The absorber surface was formed by coating the copper plate with a 
carbon black-silicon dioxide paint. The collector window was 1/8- 
inch thick glass placed one-half inch away from the collector plate. 

A 3 -mil chromel-alumel thermocouple was placed near the outer sur- 
face of the glass and shielded from radiation by covering the cavity 
in which the thermocouple lay with a small dot of white paint. 
Thermocouples were also placed on the collector plate and at the inlet 
and outlet positions of the flow channel (fig. 5). A thermocouple for 
measuring the ambient temperature was placed a few feet away from the 
collector. The overall flow loop used for the collector tests is 
shown in figure 6. A gear pump driven by a variable speed drive was 
used to control the flow of water through the collector and a rota- 
meter used for flow determination. Water at a controlled temperature 
was provided by having a temperature controller vary the power to an 
electric heater in the water tank. The radiant energy delivered by 
the simulator was measured with a water -cooled pyrheliometer. The 
emf signals from the pyrheliometer and thermocouples were recorded on 
an electronic self-balancing recorder. 

The solar simulator composed of twelve 300-2att lamps was placed 
approximately ten feet from the collector. This simulator (ref. 8) 
was used for checking out the design of the large simulator described 
in the previous section. The simulator was too small to give a 
uniformity of ±5 percent on a one square foot surface, but did provide 
±10 percent uniformity of the radiant flux. 
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Procedure 


All test runs were made at a constant flow rate (G) of 33 lb/hr 
ft • with the collector in a vertical position. The simulator was 
placed in a vertical position so that the radiation incident on 
the collector would have a zero incident angle. Once a controlled 
inlet temperature was established the power to the simulator was 
varied to determine collector performance at different radiant 
flux levels. For each adjustment of the simulator power level, 
a time of ten minutes was found to be sufficient to insure thermal 
equilibrium. Once equilibrium was determined from a visual inspec- 
tion of -the electronic pen recorder, a record of the following 
quantities was made 2 



5. Ti -T r- (water-ice mixture used for reference conditions) 

J- rex. 

6 ‘ ®dr 

The radiant flux was varied from 117 to 299 BTU/hr ft and the inlet 
temperature varied from 73°F to 109 °F. The collector efficiency was 
calculated with the following equation: 


*1 = G C p (T 0 - T-^/q^ 


( 10 ) 


Results and Discussion 


The collector performance data taken under simulated conditions 
was plotted according to equations 1, 2, and 3 as shown in figures 
7a,. 7b and 7c. The correlations are good, especially when one con- 
siders the small size of the collector and the marginal radiation 
uniformity produced by the simulator. The slope of the correlating 
line of figure 7a represents the overall heat loss coefficient (Up, ) 
which is constant for the range of absorber plates temperatures en- 
countered. A constant heat loss coefficient for the range of tem- 
peratures of the present experiments is consistent with the collector 
heat loss results of Ostad-Hosseini (ref. 9) . When radiation losses 
become more evident at higher plate temperatures than those of the 
present experiments, the overall heat loss coefficient should be- 
gin to exhibit non-linearity. The heat loss results of ref. 9 
show this to occur. 

According to equations 1, 2, and 3 the intercepts and slopes 
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of figures 7a, 7b and 7c represent performance parameters. These 
parameters are tabulated in Table I and compared to the analytical 
results of Bliss (ref. 10) for the simple collector tested. The 
comparison between the experimental and analytical results of 
Table I is good. It is clear that an analytical approach is very 
successful for a simple collector system. Difficulty in evaluating 
collector performance is presented with more complex collector de- 
signs where a more sophisticated analysis is required. For more 
complex collector designs, tests as outlined in this paper will 
conveniently and accurately give 1h e performance parameters. Gupta 
and Garg (ref. 5) have shown that the method of obtaining perform- 
ance parameters by the correlating technique of figure 7 can be 
used for most types of solar collectors. 

Some comments on the collector design may be made based on 
the performance parameters listed in Table I. The high values of 
the plate efficiency factor (F T ) and plate heat-removal efficiency 
( ) indicate a good design for efficient heat collection by the 

heat transfer medium (water in this case) . It would appear that 
there is room for improvement in the transmission qualities of the 
glass. A very clear glass can give a value of the product of ab- 
sorptivity and transmittance ( OCT e ) of about .86. The factor 
which makes this collector design poor is the overall heat loss 
coefficient ( ) which is unacceptably high. 

The overall heat loss coefficient of Table I is higher than 
one would expect from a calculation of convective losses, radia- 
tion losses, edge losses and rear losses. The higher loss en- 
countered is possibly due to the copper headers on each side of 
the collector plate. If we assume that all the collector heat 
losses (radiation and convection) are via the collector window, 
it is possible to make an estimate of the wind coefficient (h 0 ). 
This assumption is made on the basis that the collector is well 
insulated. Total heat losses are calculated by the use of equa- 
tion 1, where (Tp-T a ) . Therefore, the total heat loss equa- 

tion is as follows: 


q L = q dr (aT e C11) 

The total heat loss at the outside surface of the glass window 
is composed of convective teat losses and radiation losses as 
follows : 


1 4 f Li 

q r = h (T -T ) + e^cOT - T ) 
L o v g a J g g a J 


( 12 ) 


The radiation losses were calculated to be about 10-15% of the total 
and were subtracted from the total heat loss to give the convective 
heat loss ( q L ). The convective heat loss is plotted versus the 
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glass temperature (assuming a constant ambient temperature) in 
figure 8. The slope of the correlating line gives a wind co- 
efficient of M.O BTU/hr ft^ °F. This coefficient is equivalent 
to a wind of ten miles an hour which appears to be reasonable. 

The use of the wind coefficient has more significance than a 
statement of air speed, because of the difficulty in simulating 
wind conditions in indoor tests such as these. 

The wind coefficient and the conditions of flow, ambient tem- 
perature and incident angle are all required for correctly stating 
collector performance and for comparing the performance of differ- 
ent collector designs. 


CONCLUSIONS 

This paper has endeavored to demonstrate the advantage of 
indoor testing of flat-plate collectors under simulated conditions. 
It is believed that this type of testing allows: 

1. Convenient and standard conditions of radiation, ambient 
temperature, flow rates and wind speeds so as to permit a defini- 
tive evaluation and comparison of the performance of different 
collector designs. 

2. A determination of the key collector parameters that can 
be used for the comparative evaluation of collector designs and 
the prediction of the outdoor performance of collectors. 
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SYMBOLS 


Heat capacity, BTU/lb °F 

Collector plate efficiency factor, dimensionless 
Collector plate heat-removal efficiency, dimensionless 
Flow rate of collector fluid, lb/hr-sq ft of collector 
Free convection coefficient, BTU/hr-ft 2 °F 
Fluid convection coefficient, BTU/hr-ft 2 °F 
Wind coefficient, BTU/hr-ft 2 °F 

2 

Incident diffuse solar radiation, BTU/hr-ft 

2 

Incident direct solar radiation, BTU/hr-ft 

Incident total solar radiation, BTU/hr-ft 2 

Collector heat loss, BTU/hr-ft 2 

Collector fluid temperature, °F 

Fluid inlet temperature, °F 

Fluid outlet temperature, °F 

Collector plate temperature, °F 

Ambient temperature, °F 

Glass temperature, °F 

Effective sky temperature, °F 

Overall collector heat loss coefficient, BTU/hr-ft 2 °F 

Collector surface absorptivity, dimensionless 

Emissivity of collector surface 

Emissivity of glass 

Collector efficiency, dimensionless 



0 i Angle of incidence 

T e Transmittance 

o Stefan-Boltzmann constant, BTU/hr-ft^ °R^ 


Superscripts : 

Average quantity 

1 Absolute temperature, °R 
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Figure 4. - Air flow system for flat plate collector tests. 
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COLLECTOR EFFICIENCY, 
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RATIO TEMP DIFFERENCE TO HEAT FLUX, (Tp - T a )/qjXl0 2 , 
(BTU/HR-FT 2 -^)' 1 


(A) TEMP DIFFERENCE BASED ON PLATE TEMP. 



(0 TEMP DIFFERENCE BASED ON INLET TEMP. 

Figure 7. - Collector efficiency as function of the ratio of tem- 
perature difference to heat flux. T a , 73° F, G, 33 Ib/hr-ft 2 ; 
h 0 , 4.0 BTUfhr-ft 2 -°F. 
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